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\Q \ Abstract. The semi-relativistic quark potential model is surprisingly powerful for heavy-light 

systems if the bound state equation is treated correctly using 1 jmq expansion with heavy quark 
mass mq. We elucidate the reasons why our semi-relativistic model succeeds in predicting and 
' (— I 'I reproducing all the mass spectra of heavy-light systems so far reported, D/DjB/Bs, by reviewing 

pi I and comparing recent experimental data with the results of our model and others. Especially the 

I ' mass spectra of the so-called D^j, i.e., D*q and D' j, are successfully reproduced only by our model 

but not by other models. 
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Beginning from the discovery of the narrow meson states D*q(2317) and D'^^ (2457) (the 
so-called D^j) by BaBar [1] and CLEO [2] in 2003, respectively, open charm/bottom 
(3 ■ hadrons of the heavy-light systems have been discovered one after another. Ten years 

OO . before this discovery, we proposed a formulation for the semi-relativistic potential model 

dll], based on which we have calculated mass spectra of higher states of the heavy- 
light mesons. Subsequently to the discovery of Dsj, another set of broad heavy mesons, 
^ . Z)q(2308) and Dj (2427), were discovered by the Belle collaboration [3]. These mesons 

are identified as cq (q = u/d) excited = 1) bound states and have the same quantum 
numbers, = 0+ and 1+, as Dsj, respectively. The decay widths of these excited Dsj 
mesons are narrow, since the masses are below the DK/D*K threshold, and hence the 
dominant decay modes violate the isospin invariance, whereas the excited D mesons, 
Z)q(2308) and Dj(2427), are broad because there is no such restriction as in D^j cases. 
More recent experiments reported by CDF and DO fT] found narrow B and Bs states with 
£=1, 5i(5720), 5^(5745), (5839), and 5,i(5829). These are narrow because they 



Talk presented by T. Matsuki at the Workshop on "Scalar Mesons and Related Topics (Scadron 70)" 
held at Instituto Superior Tecnico (1ST), Lisbon, Portugal, February 1 1-16, 2008. 
^ E-mail: matsuki@tokyo-kasei.ac.jp 
^ E-mail: morii@kobe-u.ac.jp 

E-mail: kazutaka.sudoh@kek.jp 



decay through the D-waves. 

Utilizing our semi-relativistic potential model, we have so far solved the following 
problems. 

1 . Construct the formulation how to calculate the mass spectra of heavy-light systems. 

2. Numerically calculate the mass spectra of these systems and compare them with 
the experiments. 

3. Predicted mass spectra for D;q(2317), D^i(2457), D*(2308), and D\(2427) [H] 
agrees well with the experiments. 

4. Predict that 0+ and 1+ of B, are also below the threshold BK/B*K. fd\ 

5. Refurbish the calculations of [5] and fit these with the experimental data, Bi (5720), 
5* (5745), and 5^2(5839) together with the above data. |0] This calculation pre- 
dicted M{B'^^) = 5831 MeV while the experiment observes it at 5829 MeV. [4] 

6. Fit our calculations with the experimentally observed radial excitations, n ~ 2 
D;(2715) and D^(2860), and to obtain other radial excitations of D/D jB/B,. M\ 

7. Explain the superficially recovered global SU (3) invariance among 0+ states of D 
andD,. [9] 

8. Calculate the KM matrix elements by first calculating the Isgur-Wise functions 
from the wave functions used in computing the above mass spectra. [llOn 

Note that the difference between the experimental data of D^j and the threshold 
DK/D*K is only about 30 MeV but that the difference between our calculations for Bsj 
and the threshold BK/B*K is about 200 MeV. Hence the trials to explain D^j as a loosely 
bound D and K molecule can not be applied to the case for Bsj. Because it is hard to 
imagine that Dsj and Bsj have different structures, we believe that our explanation for 
these states as Qq states is legitimate both for D^j and B^j. 



OUR SEMI-RELATIVISTIC POTENTIAL MODEL 

Our formulation [5] using the Cornell potential is to expand Hamiltonian, energy, and 
wave function in terms of 1 / mg and sets coupled equations order by order. The non- 
trivial differential equation is obtained in the zeroth order, which gives orthogonal set 
of eigenfunctions, and quantum mechanical perturbative corrections to energy and wave 
functions in higher orders are formulated. Applying the Foldy-Wouthuysen-Tani (FWT) 
transformation to a heavy quark and the Hamiltonian, eigenvalue equation becomes 

H\j/e = E^\j/e, 

(//_i+//o+^i + ---)(V^£o + V^£i + ■■•)= (4' + ^i + ---)(V^^o + Wi + ---),(l) 
with the Cornell potential given by 

where integers of subscripts and superscripts denote order in 1 / mg and H = i/pwT — 
dl]. The FWT transformation is not a simple non-relativistic reduction but it also in- 
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FIGURE 1. Procedure how the degeneracy is resolved in our model. 



eludes the effects of the negative components of the heavy quark. We have the following 
expanded Hamiltonians: 

-(l+/3e)m2, (2) 
aq-p + [5q{mq + S)+V, (3) 

Here superscripts ++ mean that the matrix elements of the Hamiltonian are taken 
between the positive energy components of the heavy quarks. Negative components of 
the heavy quark have not much contributions to the masses and so are the second orders 
in 1/mg. Equation (2) gives the projection operator and determines the lowest order 
wave function which has only the positive component of the heavy quark while the 
light-antiquark is treated as a fully relativistic Dirac particle, being expressed by Eq. (3). 
The original Hamiltonian has the heavy quark symmetry (HQS) in the limit of niQ oo, 
and then this symmetry is broken by including the 1/mg correction terms. Actually the 
HQS is broken by the third term in Eq. (4), which only depends on the quantum number 
k that determines whether the HQS is broken or not. This term includes the Dirac matrix 
OCq, which has only off-diagonal matrix elements so that there is no counter term after the 
non-relativistic reduction. The chiral symmetry is broken in the first step (1), which is 
included in the Hamiltonian in a certain limit as shown in Fig. 1, then the system breaks 
the HQS in the last step (2) in Fig. 1, which is nothing but the hyperfine splitting owing 
to the third term of Eq. (4). The dominant term for the mass is given by the recoil term, 

p^/ {IniQ), the first term in Eq. (4). 
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Mass Spectra 

To demonstrate that how good our model calculations are, let us show in the fol- 
lowing Tables [2H5] the comparison with the experimental data with the parameter set 
given in Table [T] In Tables [2] ~ [51 7^ stands for the total spin and parity, Mq the low- 
est degenerate mass, c\/Mq the first order correction, Mcaic calculated value of mass, 
and Mobs observed mass. The calculated masses, M^aic, are within one percent of accu- 
racy compared with the observed masses. Mobs, k the quantum number of the operator 



TABLE 1. Optimal values of parameters. 



Parameters 

0.261 ±0.001 
mu,d (GeV) 
0.0112±0.0019 


a'' 

0.393±0.003 
m, (GeV) 
0.0929±0.0021 


a (GeV-') 
1.939±0.002 

mc (GeV) 
1.032±0.005 


b (GeV) 
0.0749±0.0020 
nih (GeV) 
4.639±0.005 





TABLE 2. D meson mass spectra (units are in MeV). 



State (2'+lLj) k Mq Ci/Mq Mcalc Mobs 



^So -1 0- 1784 0.476 x10-1 1869 1867 

^Si -1 1- 1.271x10-1 2011 2008 

1 0+ 2067 1.046 x10-' 2283 2308 

"M" 1 1+ 1.713x10-' 2421 2427 

"'Pi" -2 1+ 2125 1.415x10-' 2425 2420 

^2 -2 2+ 1.618x10-' 2468 2460 

^Di 2 1- 2322 1.894 x10-' 2762 

"^D2" 2 2- 2.054 x10-' 2800 



—Pq \ Lq • L + 1 j , which denotes the degenerate states. In the calculations, we have used 
values of parameters listed in Table 1 . 

Comparison with Other Models 

After observing that our model nicely succeeds in predicting and/or reproducing the 
experimental data for the heavy-light mesons, we should clarify the reason why our 
model well works while others do not. Especially the other models have trouble to 
generate masses for the 0+ and 1+ states of D^. We will give Table|6] which qualitatively 
describes the differences between our model and others. As one can see in Table [6l 
only successful quark potential model to reproduce masses of D^j is our semirelativistic 
model. A coupled channel method is also successful but the physical meaning remains 
obscure in that the authors of [|l5il do not take into account all the channels. 

The BS equation is proposed by Zeng, Van Orden, and Roberts [12] to describe 
the heavy-light system, which is similar to ours except that they neglect the negative 
components of the heavy quark. Their numerical calculations give values higher than 
DK/D*K thresholds and use constituent quark masses. The differences between ours 
and theirs are i) whether the light quark masses niq are small or not, i.e., current or 
constituent quark masses, ii) whether the negative components of the heavy quark are 
taken into account or not. We adopt the current quark masses, mj, = = 11.2 and = 
92.9 MeV while they adopt iriu = nid = 248 and = 400 MeV. We take into account 
the negative components of the heavy quark which contribute to the second order 
calculations in 1/mg, while the paper [12] takes into account the second orders coming 
from only the positive components of the heavy quark. Considering our successful 



TABLE 3. Dj meson mass spectra (units are in MeV). 



State C'+'^Lj) 


k 


7^ 


Mq 


ci/Mo 


Mcalc 


Mobs 




-1 


0- 


1900 


0.352 X 10-' 


1967 


1969 




-1 


1- 




1.102 xlO-' 


2110 


2112 




1 


0+ 


2095 


1.101 xlO-' 


2325 


2317 


"M" 


1 


1 + 




1.779 xlO-' 


2467 


2460 


"iPi" 


-2 


1 + 


2239 


1.274 X 10-' 


2525 


2535 


'Pi 


-2 


2+ 




1.467 xlO-' 


2568 


2572 




2 


1- 


2342 


2.032 X 10-' 


2817 




"3£)2" 


2 


2- 




2.196 xlO-' 


2856 





TABLE 4. B meson mass spectra (units are in MeV). 



State C'+^Lj) 


k 




Mo 


ci/Mo 


Mcalc 


Mobs 


'So 


-1 


Q- 


5277 


-0.161 xlO-2 


5270 


5279 




-1 


\- 




0.981 xlO-2 


5329 


5325 


'Po 


1 


0+ 


5570 


0.401 xlO-2 


5592 






1 


1 + 




1.412 xlO-2 


5649 




"'Pl" 


-2 


1 + 


5660 


1.069 xlO-2 


5720 


5720 


'Pi 


-2 


2+ 




1.364 xlO-2 


5737 


5745 


'D, 


2 


1- 


5736 


2.203 xlO-' 


6999 




"3^2" 


2 


2- 




1.430 xlO-' 


6556 





calculations, we believe that if they [fl2|] adopt the current quark masses, then they would 
obtain the correct mass values for D^y by adjusting parameters. How the light quark mass 
affects the spectra can be seen in Figure 1 of Ref. [17], in which paper the average D 
meson mass of D and D* is calculated by varying the c quark mass and by taking two 
values of the light quark mass, mj, = 10 and 336 MeV. Even though the potential form is 
different from ours, this figure shows that the value of the light quark mass is important 
to determine the spectra of the heavy-light system. It turns out that only the case of 
niu = 10 MeV, i.e., current quark mass, can fit with the experiments for the heavy-light 
system. 

In the last two rows of Table [6l we list the experiments and the sum of masses of 
D/D* and K so that one can see how far calculated ones are away from the experiments 
and D/D*K thresholds. 
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TABLE 5. Bs meson mass spectra (units are in MeV). 



State (^^+'Ly) 


k 


jP 


Mn 


n / Mn 


-"^^calc 






-1 


0- 


5394 


-0.302 xlO-2 


5378 


5369 




-1 


l- 




0.853 xlO-2 


5440 






1 


0+ 


5598 


0.350 xlO-2 


5617 






1 


1+ 




1.498 xlO-2 


5682 






-2 


1+ 


5775 


0.978 xlO-2 


5831 


5829 


'P2 


-2 


2+ 




1.263 xlO-2 


5847 


5839 




2 


1- 


5875 


2.949 xlO-2 


6048 






2 


2- 




0.564 xlO-2 


5908 





TABLE 6. Comparison ours with other models. QPM means quark potential model 



Method 


Authors 


Dsj 


Successful? 




Semirelativistic QPM 


T. M. etal.[5, 7] 


2.339, 2.487 GeV 


OK 


current 


Conventional QPM 


Godfrey et al.[llj 


2.48, 2.55 


No 


constituent 


BS eq. ^ ours 


Zeng et al.[121 


2.38,2.51 


No 


constituent 


Another QPM 


Ebertet al.[13J 


2.463,2.535 


No 


constituent 


DK Molecule 


Barnes et al.[14j 




? 


N/A 


Coupled Channel 


Beveren et al.ri51 


2.28 (2.320) 


OK 


N/A 


tetraquark 


Cheng et al.[16J 




? 


N/A 


Observed 


[1,2] 


2.317,2.460 






D + K/D''+K 




2.367, 2.505 
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